An alternative normal state c-axis resistivity model for high-T c superconductors 
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An alternative model for c-axis resistivity in layered high-T c crystalline superconductors is pro- 
posed and has been characterized as an essentially two-dimensional Fermi liquid. Average ionization 
energy is included as additional parameter that determines the concentration of tunnelling electrons 
between CU-O2 layers. This model agrees well quantitatively with the Bi2212 and Y123 single 
crystals, and qualitatively with the pure 1212 phase polycrystals. 
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1. Introduction 

The out-of-plane conduction mechanism in the nor- 
mal state of high-T c superconductors is intriguing. In 
the early stages, facts were pointed out how anoma- 
lous various normal state properties were in the simple 
framework of Fermi liquid This Fermi liquid pre- 
dicted large Fermi surface whereas small Fermi surface 
was expected from the strong correlation models where 
the normal state transport properties favored the latter. 
Further studies revealed that the conventional field the- 
oretic many-body treatment does not seem to lead to 
the non- Fermi liquid behavior either Q. Apart from 
that, according to RVB theory, the normal state of high- 
Tc superconductors is all metals in the CU-O2 planes 
and is semiconductor-like for conduction between those 
planes Q. In addition, Anderson and Zou suggested that 
both the a-b plane and c-axis resistivities can be fitted 
very accurately by p = A/T + BT. It is also believed 
that there is a gap induced by the localization of carriers 
along c direction, which acts as an anomalous constraint 
of c-axis hopping rate Parallel to this, the c-axis 
normal state resistivity of high-T c superconductors were 
derived via Fermi-Dirac statistics assuming the only ob- 
ject that were allowed to tunnel through the CU-O2 lay- 
ers are the real electrons as required by the RVB theory. 
This phenomcnological model is also used to fit the a-b 
plane resistivity. Furthermore, the estimation of resis- 
tivity transition such as from metallic to semiconductor 
or vice versa upon substitution in the pure 1212 phase 
polycrstals is also reported. 



2. The model 

The tunnelling electron's distribution can be derived 
using Fermi-Dirac statistics and ionization energy as an 
anomalous constraint. This derivation involves two re- 
strictive conditions: i) the total number of electrons in 
a given system is constant and ii) the total energy of n 
electrons in that system is also constant. Both conditions 
are as given below 
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In condition (2) however, the ionization energy, Ej will 
be included as an additional constraint where E e i ectron = 
Einitiai state + £/• This is to justify that an electron 
to occupy a higher state N from the initial state M 
is more probable than from state L if the condition 
Et(M) < Er{L) at certain temperature, T is satisfied. 
Thus, condition (2) is rewritten as 
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The importance of this inclusion is that it can be inter- 
preted as a gap and also, particularly the Ej can be used 
to estimate the resistivity transition upon substitution of 
different valence state ions. Finally, the probability of an 
electron to occupy an energy level E is given by 
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if e^+KE+E,) ^ 1 Note that the distribution (4) can 
be reduced to the standard Fermi-Dirac distribution if 
Ei — > 0. The concentration of electrons is given by 



f e (E)N e {E)dE, 
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where fj, = -E F /kT, A = 1/kT and N e (E) is the two 
dimensional density of states that was derived from the 
schrodinger equation. The solution of integral (5) gives 
the concentration of tunnelling electrons between CU-O2 
layers for high-T c superconductor as 
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This concentration is expressed in term of ionization 
energy and k is the Boltzmann constant, m* is the ef- 
fective electron mass and Ep is the Fermi level. It is 
well known that the scattering among Fermions gives 
1/t oc T 2 and p oc T 2 from the empirical rule @, H, 0] • 
According to the elementary metal theory, the resistivity 
of metal is given by p = m*/ne 2 T where e and r are the 
electron's charge and the mean scattering free time, re- 
spectively Therefore, the phenomenological model of 
c-axis normal state resistivity can be written in the form 
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Where A is independent of temperature. This model is 
comparable with the form given in Ref. Q where in the 
exponential term, the gap, S is used instead of the term 
(Ei — Ep) as an anomalous constraint of c-axis hopping 
rate between Cu-02 layers. It shows that this exponential 
term has a similar meaning to this gap. 



3. Discussion 



temperature at T < 155. 7K and T < 139. 8K respec- 
tively suggesting the semiconductor-like region starts be- 
low those temperatures. The calculated p a b from model 
(7) for samples B and C as shown in Fig. 2b were also in 
good agreement with the experimental data. However, 
the Ei — Ep parameter has a much lower value as ex- 
pected compared to c-axis which can be due to the low 
spin fluctuations in the Cu-O— 2 planes. There is no sig- 
nificant differences in A and Ei — Ep parameters between 
sample B and C for both the a — b plane and c-axis al- 
though the ratio a/c differs in large magnitude, which 
are approximately 8 for B and 16 for C. 

The resistivity transition from sample la to lb as the 
O2 content decreases in Bi2Sr2CaCu20j, is due to the 
transition of Ei - E F (110K) < kT r to Ei - Ep (390K) 
> kT r , T r is the room temperature. Therefore, the con- 
centration of electrons, n between layers in sample lb 
is small due to insufficient thermal energy since kT r < 
390K. Table 1 lists A and Ei — Ep parameters in detail 
for the calculated curves. In the p(T) transition estima- 
tion, Tl-Sr-(Yi_/X/)-Cu-0 system is considered for an 
example where the average Ei for X z+ and ions in 
that system can be calculated using 



The experimental p c plots of doping dependent Bi2212 
and Y123 single crystals were best fitted quantitatively 
with the phenomenological model (7). The Bi2212 sam- 
ples, la and lb were annealed at 10 and 1CT 3 atm respec- 
tively to control the oxygen content Q. In contrast, all 
the Y123 samples have the same composition, which is 
YBa 2 Cu307, except for the volume of the unit cell which 
are 630 x 640 x 75 and 390 x 400 x 25 ^m 3 for samples 
B and C respectively @. 

The experimental plots and the calculated curves of 
p c (T,Ei) versus T for Bi2212 are shown in Fig. la and 
b. Both the experimental plot and the calculated curve 
in Fig. la can be separated into two regions: (i) Ei — Ep 
(110K) < kT (metallic-like behavior at T > 110K) and 
(ii) Ei — Ep (110K) > kT (semiconductor-like behavior 
at T < 110K). That is, p c (T,Ei) property changes from 
metallic-like behavior at high temperatures to semicon- 
ductor behavior at low temperatures. At high temper- 
atures, p c oc p a b, suggesting that scattering or fluctua- 
tions in the CU-O2 planes may dominate the linear term 
of p c @, HI • In addition, the thermal energy, kT is higher 
than the constraint, Ei — Ep above 110K which give 
rise to the conductivity of tunnelling electrons. At low 
temperatures however, p c is incoherent and can be inter- 
preted as an indication of a low thermal energy as T goes 
below 110K. 

Similarly, the experimental and the calculated resis- 
tivity curves for Y123 samples B and C as shown in 
Fig. 2a also indicates the same characteristics, where 
the metallic-like behavior occurs when Ei — Ep (155. 7K) 
< kT for sample B and Ej - E F (139.8K) < kT for sam- 
ple C. The resistivity for both B and C increases with 
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and 
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i = 1, 2,...z and i = 1, 2,. ..A. z + and A + is the valence 
state of X and Y respectively. If Ei[X z+ ] < Ei[Y A +] then 
the substitution of Y with X will reduce Ei — Ep param- 
eter thus, reduces p c (T) as the curves in Fig. 3 shows. 
Furthermore, the p c (T) rises rapidly with Ei — Ep for 
100K compared to 300K due to the low tunnelling elec- 
tron's concentration. This figure also reveals that the 
room temperature c-axis resistivity increases exponen- 
tially with the systematic substitution of different valence 
state ions. Table 2 summaries the estimated p(T) for 
single and double substituted polycrystallinc high-T c of 
pure 1212 phase superconductors. Substitution of Sr 2+ , 
Er 3 +, Pr 3+ and Y 3 + with Ba 2 +, Sr 2 +, Sr 2+ and Ca 2 + 
respectively reduces p(T) since the ionization energies of 
the latter is lower than the former. On the other hand, 
the substitution of Sr 2+ with Tm 2 +, Sm 2+ , Eu >2 + and 
Dy >2+ increases Ei hence, p(T), therefore the metal- 
lic to semiconductor-like transition can be expected and 
vice versa for the former. Here, p c (T) is assumed to be 
roughly p(T) due to polycrystallinity of the samples. 
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Figure 1: Temperature dependence of the c-axis resistivities 
of Bi2212 single crystals annealed at a) 10 and b) 10 -3 atm. 



Figure 2: Experimental (dotted) and calculated (solid line) 
plots for YBa2Cu2C>7 single crystals, a) c-axis and b) a — b 
plane resistivities. 



4. Conclusions 

Ionization energy has been used in the Fermi-Dirac 
statistics as an additional restrictive condition. The cal- 
culation for the tunnelling electron's concentration was 
done by employing the two-dimensional Fcrmion charac- 
teristics. The proposed model (7) agrees well quantita- 
tively with the experimental p c (T) of Bi2212 and Y123 
single crystals and also predicts the resistivity will in- 
crease with the average ionization energy. In addition, 
both the experimental data of p c (T) and p a b(T) can be 
fitted with the proposed model. Besides that, the es- 
timated p(T) transition upon substitution of ions with 
larger or smaller ionization energies also agrees well with 
the pure 1212 phase polycrystals. 
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Figure 3: The variation of (Ei — Ef) parameter with c-axis 
resistivity at 100 and 300K. 

Table I: Detailed values of temperature independent A and 
Ei — Ef parameters. 
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pure 1212 phase polycrystalline high-T c superconductors. 
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Samples 


A (ms 1 K ') 


Ei-Ev (J) 


la 


1.4 x 10 3 


110 K 


lb 


6.117 x 10 3 


390 K 


Bab 


5.707 x 10 3 


36.6 K 


C ab 


4.636 x 10 3 


38.6 K 


B c 


2.12 x 10 5 


155.7 K 


Q 


2.832 x 10 5 


139.8 K 
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